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ABSTRACT 

In this paper , we considered a nonlinear two-dimensional quadratic Volterra integral equation of the second kind with 
continuous kernel. The existence of a unique solution of this equation under certain conditions, was discussed in the 
Banach space C ([0 , T] x [0, T]) using Pickard method. Adomian method was applied on the Two-dimensional quadratic 
Volterra integral equation for a nonlinear case. Two applications were presented; Maple 18 software was used to 
compute the numerical solutions. Finally, we realized that the Adomian method is appropriate and effective for treating 
two-dimensional linear and nonlinear quadratic integral equations. Moreover , we noted that the effect of time factor was 
obvious on numerical results. 
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1. INTRODUCTION 

Quadratic integral equations (QIEs) include several types of integral equations. QIEs are often applicable to the 
theory of neutron transport, traffic theory, modelling of radiative transfer, the queuing theory, the kinetic theory of 
gases to name a few amongst many other phenomena. 

Chandrasekhar equation is an important QIE. This equation describes a radiative transfer. An important 
special case of that equation is Chandrasekhar’s integral equation (CIE) which appears in radiative transfer. There 
are three areas where CIE is applied which includes stellar atmosphere, the kinetic theory of gases and neutron 
transport, see Chandrasekhar [1,2], Argyros [3], Caballero et al., [4], Hu et al., [5], and Cahlon [6]. 

A number of researchers have an interest in QIEs. Concerns lie majorly in numerous studies on the 
analytical properties of quadratic integral solutions; the properties of QIEs include compactness, the existence of 
solutions for several classes, monotomic and positive solutions. Abbas et al., [7], Hashem and Alhejelan [8], 
Darwish [9] and Banas et al., [10] used the theory of measures of noncompactness to prove the existence of 
theorems for QIEs. Wang et al., [11] studied the existence of positive solutions for the NQIE. El-Sayed and 
Hashem [12] presented an existence theorem for at least one continuous solution for NQIEs of fractional order. 
Darwish and Ntouyas [13] presented an existence theorem for monotonic solutions. Zhu [14], using the theory of 
measures of noncompactness. 

Some numerical and analytical methods can be applied to estimate the solutions of the QIEs. However, 
ADM is the most common method used to obtain the numerical solutions for linear and nonlinear QIEs. In El- 
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Sayed et al., [15] the Picard and ADM methods were applied in QIEs. Abdou and Raad [16], Raad [17] used the modified 
ADM, with the aid of modified Simpson's rule (MSR). 

2. TWO-DIMENSIONAL VOLTERRA QUADRATIC 
INTEGRAL EQUATION OF THE SECOND KIND 

Consider the T-DVQIE of the second kind 

t ) = f(t,T)+x(t,T,<p(t,T))f*^k(t,r-,T,s)Y(r,s,<p(r,s))drds. (1) 

In (1) ju is a constant and /i ^ 0, the functions y( r,s,<p( r, s)) and (t, x, <p(t, r)) „ are known as continuous 
nonlinear functions. k(t,r;x,s) represents the kernel of the integral equation which will be continuous and the 
function/(t, r) is a known linear continuous function, while 0(t, r) is unknown, will be determined in the 
spaceC([0, 7] x [0,7]). 

3. EXISTENCE OF A UNIQUE SOLUTION OF T-DVQIE 

In this section, the existence of a unique solution of equation (1) will be discussed and proved using Picard theorem. 
Considering this, we assume the following conditions: 

• The kernel k(t, r; x, s)e C([0, 7] x [0,7])is a continuous function and satisfies: 

| k ( t , r;x,s)\ < L, (H is constant). 

• The given function / ( t , r) and its partial derivatives with respect to t, x are continuous in C([0, T] X [0, T]). 

• The nonlinear function y(t, r, $(t,r)) is continuous and satisfies for the constants A 1: A 2 the conditions: 

(a) \y(t,T,(f>(t, t))| < A v 

(h) | y(t, t, 0!(t, t)) - y(t, t, (p 2 (t, t)) |< A 2 \ fa (t, r) - (p 2 (t, r) | . 

• For the constants B 1 ,B 2 , The nonlinear function j(t, r, (p(t, r)) is continuous and satisfies the conditions: 

(a) |j(t, t, (p(t, t))| < B x . 

(b) | x(t, t, 0!(t, t)) - x(t, t,( p 2 (t,r))|< B 2 \cp 1 (t,r) - (p 2 (t, t)| . 

Theorem 

The solution of (1) exists and is unique in the Banach space C([0, T] X [0,7]), under the condition \jj \ > L [A 2 B 1 + A 1 B 2 \. 

Proof 

Using Picard method, by picking up any real continuous function 0 o (t, r) in the spaceC([0,7] X [0,7]), then 
constructing the solution as a sequence of functions {0 n (t, x )},n -> oo. 

For every 0 n , equation (1), will be in the form 

B<Pn(t, t) = fit, t) + x(t,T,(pn-i(t,T)) I* I* k(t,r-,T,s)y(r,s,(p n _ 1 (r,s))dsdr, 

n = 1, 2, 3, 
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n 0 o (t, t) = / (t, t), (2) 

functions 0 n (t, r) are continuous and can be written as a sum of differences: 

0n = 0 o + X;=i(0j - 4>j-i) (3) 

This means that the convergence of the sequence 0 n is equivalent to convergence of the infinite series and the 
solution will be 




(pit, t) = t), (V f,r £ [0,7]) (4) 

From (2) for n = 1, we get 

|0i(t, T)-0 o (t. t)| = |i| |/(t,T,^ 0 (t,T))|/ 0 t / 0 T |fc(t,r;r,s)| |y(r,s,0 o (r,s))|dsdr, 
using assumptions (i), (iii — a) and (iv — a) 

101 (t, T)- 00(t, r)| < |^| i Jo J 0 T dsdr < |^-| it t. (5) 


Now, we shall obtain an estimate for n > 1, from (2) we can get 
M[0„(t,r)-0 n _i(t,r)] = 

x(t, t, 0 n _ 1 (t, t) - 0„_ 2 (t, t)) / 0 C / 0 t fc(t, r; r, s)y(r, s, 0 n _ 1 (r, s) - 0 n _ 2 (r,s)) dsdr, 
taking the absolute value of equation (6), to obtain 

|0„(t, t)- 0 n _ 1 (t, t)| = |i| |j(t,r,0 n _ 1 (f,r))|/ o t / o T |fc(t,r;T,s)||y(r,s,0 n _ 1 (r,s))|drds 

- |^| |t, 0„_ 2 (t, t))| Jq / 0 T |fe(t, r; r, s)||y(r,s, 0 n _ 2 (r,s)) |drds. 

When n = 2, we have 
|0 2 (t,r) - 0 1 (t,r)| < |i| L 2 i4 1 B 1 f 2 r 2 

repeating this technique, we obtain the general form for the terms of the series 

l0n(t.T)-0„_i(t, t)| < A l B l [^A 2 B 1 + A^] 

x [^ 2 Bi + A ± B 2 ]... [^z^A 2 B 1 + A^] [^A 2 B 1 + A x B 2 

< (|i| L t t [A 2 B ± + Ajff < p n , (p = |i| L [A 2 B 1 + A^]). 

According to the condition ( L [A 2 B X + AjB 2 ] < |(i|), when n -»oo, the sequence {c(>](t,x)} is uniformly convergent, 
and 


( 6 ) 


(7) 


(8) 


0(t,T) = T,t =0 4>i(t, t). 


(9) 
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Since each of 0;(t, t) in (9) is continuous, therefore 0(t, r) is also continuous. 
Thus, the existence of a continuous solution of equation (1) is proved. 

To prove 0(t, t) is unique, assume cp{t, r) is another continuous solution, i.e. 
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n(p{t, r) = fit, t) + x(t,T,<p(t,T)) f* f* k(t,r;T,s)Y(r,s,(p(r,s))dsdr, (10) 

hence, we get 

\(p(t,T) - cp n (t,T)\ < |i| \xit,r,(p)\ f* Jj | k(t, r; r,s)11y(r, s, cp) - yir,s, (p n _ i)| drds 

+ |^| I xit, T, (p) - xit, T, (pn-f )| Jq J 0 t | kit, r; t, s ) 11 Yir, s, (pn^drds. 

Using conditions (iii— a) and (iv—a), with the aid of equation, we get 

\(pit,T) - (p n it,r)\ < p n , p < 1, (11) 

hence 

lim n^cfnit, r) = (pit, r) = (pit, r), (Vt,t £ [0,T]). 
which prove that the solution (pit, r) of (1) is unique. 


4. ADOMIAN DECOMPOSITION METHOD FOR TWO DIMENSIONAL 
NONLINEAR VOLTERRA QUADRATIC INTEGRAL EQUATION 


ADM has been known as a reliable scheme for solving a variety of linear and nonlinear integral equations. In many papers 
ADM has been used to solve some classes of integral equations in one dimension. The Adomian polynomials guarantee the 
convergence of the solution. The convergence of ADM is discussed and proved by different methods. 

Here, the T-DNVQIE (1) will be solved using ADM. For this, we express the solution of (1) in the series 

form by 

(pit.T) = T,tf 0 (pi (t, r), (12) 

the decomposition method identifies the nonlinear terms y(t, r, (pit, r)), y(t, r, (pit, r)) by 


y(t, T, 0(t, t)) = K=0^n- x(t,T,(p(t, t )) = Y.m=0 A m , 

where the Adomian polynomials A n , A m can be determined by the relation: 


p =o 


p=o 


Application of ADM to equation. (1) yields 

/j(poit,r) = fit, t ); 


(13) 


(14) 


A Kpiit, t) = fit, t) + A m fg fjkit, r; r, s') A n drds, (t > 1). 


(15) 


5. NUMERICAL RESULTS 
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Volterra Integral Equation of the Second Kind 

Consider the T-DQVIE (1), when p = 1 „ in three different times. The time interval will be divided by n = m = 10. 
The values in the tables are mentioned only in t = x. The following tables give an exact solution, the approximate solution 
(Num. Sol.), and their corresponding errors (Error) of both linear and nonlinear types of y(r, s, 0(r, s)). The diagrams 
explain the difference between these results. 

Application (1) 

Consider the T-DQVIE 

(pit, T) = fit, t) + ^ W> 3 (t,r)] jgfg(2r + t + s - r) (p k {s,r ) drds. (16) 

This application was solved for the linear type when k = 1, and k = 2 for the nonlinear type. We took three-time 
values: T = 0.067, 0.1, 0.45. 

The exact solution is 0(t, t) = sin (t + r)) 


Table 1: t = x = 0.067 


T 

T 

Exact 

Linear Nonlinear 

Num. Sol. 

Error 

Num. Sol. 

Error 

0 

0 

0 

0 

0 

0 

0 

0.0134 

0.0134 

0.042084909 

0.042084908 

0 

0.042084909 

0 

0.0268 

0.0268 

0.084095246 

0.084095246 

6.000E-11 

0.084095246 

1.000E-11 

0.402 

0.402 

0.125956573 

0.12595657 

1.000E-09 

0.125956573 

3.000E-10 

0.0536 

0.0536 

0.167594714 

0.16759472 

9.500E-09 

0.167594717 

2.400E-09 

0.067 

0.067 

0.20893589 

0.208935934 

4.460E-08 

0.208935904 

1.400E-08 


Table 2: t = x = 0.1 


t 

T 

Exact 

Linear 

Nonlinear 

Num. Sol. 

Error 

Num. Sol. 

Error 

0 

0 

0 

0 

0 

0 

0 

0.02 

0.02 

0.06279052 

0.06279052 

1.000E-11 

0.06279052 

0 

0.04 

0.04 

0.125333234 

0.12533323 

1.300E-09 

0.125333234 

3.000E-10 

0.06 

0.06 

0.187381315 

0.18738133 

2.070E-08 

0.187381321 

5.900E-09 

0.08 

0.08 

0.248689887 

0.24869004 

1.518E-07 

0.248689944 

5.680E-08 

0.1 

0.1 

0.309016994 

0.30901770 

7.031E-07 

0.309017321 

3.269E-07 


Table 3: t = x = 0.45 


T 

T 

Exact 

Linear Nonlinear 

Num. Sol. 

Error 

Num. Sol. 

Error 

0 

0 

0 

0 

0 

0 

0 

0.09 

0.09 

0.278991106 

0.278991447 

3.4130E-07 

0.278991249 

1.4310E-07 

0.18 

0.18 

0.535826795 

0.535862705 

3.5910E-05 

0.535855961 

2.9166E-05 

0.27 

0.27 

0.75011107 

0.750544338 

4.3327E-04 

0.750612699 

5.0163E-04 

0.36 

0.36 

0.904827052 

0.906705029 

1.8780E-03 

0.907550904 

2.7239E-03 

0.45 

0.45 

0.987688341 

0.991512576 

3.8242E-03 

0.994225947 

6.5376E-03 
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t = t = 0.067 






The previous numerical results of Tables (1) to (3), have shown that 


• It is clear the convergence of the approximate solution to the exact solution. 

• The correlation between the error and the variables t and r. 

• Nonlinear type results are better than linear type in small time periods. 


Application (2) 

Consider the T-DQVIE 

0(t,r) = /(t,r) + ln(0(t,r)) x /JCrt + ST)y(r,s,(p(r,s )) drds, (17) 

where in linear type y(r, s, 0(r, 5 )) = 0(5, r), and the nonlinear type y(r, s, 0(r, 5 )) = exp(0(5, r)), we pick three-time 
values: T = 0.009, 0.07, 0.25 and the exact solution is 0(t, r) = (t + r + 1) 


Table 4: t = x = 0.009 


T 

T 

Exact 

Linear 

Nonlinear 

Num. Sol. 

Error 

Num. Sol. 

Error 

0 

0 

1 

1 

0 

1 

0 

0.0018 

0.0018 

1.00360000 

1.00360000 

0 

1.00360000 

0 

0.0036 

0.0036 

1.00720000 

1.00720000 

0 

1.00720000 

0 

0.0054 

0.0054 

1.01080000 

1.01080000 

0 

1.01080000 

0 

0.0072 

0.0072 

1.01440000 

1.01440000 

0 

1.01440000 

0 

0.009 

0.009 

1.01800000 

1.01800000 

0 

1.01800000 

0 


Table 5: t = x = 0.07 


T 

T 

Exact 

Linear 

Nonlinear 

Num. Sol. 

Error 

Num. Sol. 

Error 

0 

0 

1 

1 

0 

1 

0 

0.014 

0.014 

1.02800 

1.028000000 

0 

1.028000000 

0 

0.028 

0.028 

1.05600 

1.056000000 

0 

1.056000002 

2.0000E-09 

0.042 

0.042 

1.08400 

1.084000009 

9.0000E-09 

1.084000025 

2.5000E-08 

0.056 

0.056 

1.11200 

1.112000049 

4.9000E-08 

1.112000144 

1.4400E-07 

0.07 

0.07 

1.14000 

1.140000183 

1.8300E-07 

1.140000553 

5.5300E-07 
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Table 6: t = x = 0.25 


T 

T 

Exact 

Linear 

Nonlinear 

Num. Sol. 

Error 

Num. Sol. 

Error 

0 

0 

1 

1 

0 

1 

0 

0.5 

0.5 

1.1 

1.100000025 

2.50000E-08 

1.100000073 

7.30000E-08 

0.1 

0.1 

1.2 

1.200001517 

1.51700E-06 

1.200004782 

4.78200E-06 

0.15 

0.15 

1.3 

1.300016503 

1.65030E-05 

1.300055771 

5.57710E-05 

0.2 

0.2 

1.4 

1.400088306 

8.83060E-05 

1.400317411 

3.1741 IE-04 

0.25 

0.25 

1.5 

1.500318742 

3.18742E-04 

1.501201723 

1.20172E-03 








The previous numerical results of Tables (4) to (6) have shown that: 

• The Adomian method can be flexible to effectively solve various problems. 

• The effect of the time factor is evident on the numerical solution. 

• The results of the linear type are better than the nonlinear type. 

6. CONCLUSIONS 

By the current research, we considered a two-dimensional quadratic Volterra integral equation of the second kind with 
continuous kernel. The successive approximation method has been used to prove the existence and uniqueness solution. 
Moreover, two applications have been solved to illustrate the efficiency of the ADM. The effective effect of ADM 
application on quadratic integral equations is obvious. 

Several cases can be established when the T-DVQIE takes special forms. For example, the T-DVQIE of the first 
kind and when the kernel k(t, r; t, s ) is discontinuous in the future work. 
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